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ABSTRACT 

Cosmological fluids are commonly assumed to be distributed in a spatially homogeneous 
way, while their internal properties are described by a perfect fluid. As such, they influence 
the Hubble-expansion through their respective densities and equation of state parameters. The 
subject of this paper is an investigation of the fluid-mechanical properties of a dark energy 
fluid, which is characterised by its sound speed and its viscosity apart from its equation of 
state. In particular, we compute the predicted spectra for the integrated Sachs-Wolfe effect 
for our generalised fluid, and compare them with the corresponding predictions for weak 
gravitational lensing and galaxy clustering, which had been computed in previous work. We 
perform statistical forecasts and show that the integrated Sachs-Wolfe signal obtained by cross 
correlating Euclid galaxies with Planck temperatures, when joined to galaxy clustering and 
weak lensing observations, yields a percent sensitivity on the dark energy sound speed and 
viscosity. We prove that the iSW effect provides strong degeneracy breaking for low sound 
speeds and large differences between the sound speed and viscosity parameters. 

Key words: cosmology: weak gravitational lensing, cosmic microwave background, meth¬ 
ods: analytical 


1 INTRODUCTION 

The expansion dynamics of the Universe is usually described by as¬ 
suming ( i ) general relativity as the theory of gravity, (if) a high de¬ 
gree of symmetry, namely spatial isotropy and homogeneity at each 
instant in time, and (Hi) ideal fluids which source the gravitational 
fields. These three assumptions lead to the Friedmann-equations for 
the time-evolution of the scale factor a(t), which reflect the fact that 
Einstein’s field equation is of second order, and shows acceleration 
or deceleration a as phenomena. 

The inclusion of the cosmological constant on grounds of 
the Lovelock theorem, which states that the field equation is the 
most general one in four dimensions, which includes derivatives 
of the metric of up to second order, and which conserves energy- 
momentum, yields a natural way to explain cosmic acceleration at 
late times. 

Introducing dynamic dark energy components based on scalar 
self-interacting fields and interpreting the energy-momentum- 
tensor with the corresponding conservation law allows the iden¬ 
tification of the homogeneous and isotropic field with a relativis¬ 
tic ideal fluid, whose relation between pressure and density is 
parametrized by an equation of state w = p/p. At the same time, 
this equation of state is the only free function that is allowed by 
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the Einstein field equation with the symmetry assumptions of the 
Robertson-Walker metric. 

For these reasons, a central goal of cosmology is to investigate 
dark energy and the cosmological constant through their influence 
on the dynamics of the scale factor and on the growth of structures. 
The fluid-picture is attractive due to its generality: Apart from ac¬ 
tual substances like relativistic components (w = +1/3) and nonrel- 
ativistic components (w = 0), it is general enough to describe spa¬ 
tial curvature (w = -1/3) and the cosmological constant (w = -1), 
while isotropy and homogeneity of the fluid ensure the Friedmann- 
symmetries. 

Dark energy models based on self-interacting scalar fields 
show a natural variation of the dark energy equation of state pa¬ 
rameter, because their time-evolution is governed by the Klein- 
Gordon equation, and therefore the kinetic and potential terms in 
their energy-momentum tensor evolve, leading to a time evolution 
in the equation of state, and therefore to a variation of their influ¬ 
ence on the expansion dynamics of the Universe. In the slow-roll 
limit one recovers values of w close to -1, resulting in accelerated 
expansion. 

Adopting the fluid picture is an important test of whether the 
dark energy fluid is ideal or not: If the fluid has inhomogeneities, 
pressure fluctuations and density fluctuations are related through a 
sound speed, which in the most straightforward case describes an 
adiabatic compression of the fluid, there can be anisotropic stresses. 
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and finally, velocity perturbations can experience viscous forces 
that dissipate kinetic energy (for literature in this field we refer 
to Battye & Moss 2006; Mota et al. 2007; Battye & Moss 2009; 
Calabrese et al. 2011; Ballesteros et al. 2012; Sapone & Majerotto 
2012; Sapone et al. 2013; Appleby et al. 2013; Dossett & Ishak 
2013; Sawicki et al. 2013; Amendola et al. 2014; Chang & Xu 
2014; Chang et al. 2014; Cardona et al. 2014; Pearson 2014; Balles¬ 
teros 2015). In addition, there can be a nonlinear relation between 
pressure and density of a fluid, one example of which would be 
Chaplygin-cosmologies (Bento et al. 2002; Li & Xu 2014), while a 
similar phenomenology could in principle be due to modifications 
in gravity rather than due to non-ideal fluids under general relativity 
(Kunz & Sapone 2007; Bertschinger & Zukin 2008; Silvestri 2009; 
Pogosian et al. 2010; Song et al. 2010; Leon & Saridakis 2011; 
Saltas & Kunz 2011; Baker et al. 2013; Boubekeur et al. 2014). 

These modifications break homogeneity on small scales, and 
require corresponding fluid equations for their time evolution, as 
well as couplings to local gravitational fields, which enable interac¬ 
tion between the dark energy fluid, the dark matter and the baryonic 
component. Commonly, one observes a difference between the two 
metric potentials in the case of nonzero sound speeds and equations 
of state unequal to -1, which can be probed by photons, relative to 
the motion of nonrelativistic objects such as galaxies, which is only 
sensitive to a single metric potential. 

In this paper we investigate cosmological perturbations with 
a non-ideal dark energy fluid and aim to forecast constraints on its 
speed of sound c s and its viscosity from Euclid 1 (Laureijs 2009; 
Laureijs et al. 201 1) and Planck (Ade et al. 2014b, 2015). Specifi¬ 
cally, we consider tomographic weak gravitational lensing (Ayaita 
et al. 2012), galaxy clustering (DeDeo et al. 2003; Takada 2006) 
and the integrated Sachs-Wolfe effect (Dossett & Ishak 2013; So- 
ergel et al. 2015) as probes on the influence of non-ideal fluids on 
the statistics and the evolution of structures, while the background 
expansion dynamics is given through the individual density param¬ 
eters and the equation of state parameters, assuming that there is no 
energy exchange between the fluids. 

Our work is complementary to that of Mota et al (2007); Cal¬ 
abrese et al. (2011); Chang & Xu (2014), who used the same model 
to describe the evolution of anisotropic stress: Mota et a). (2007); 
Chang & Xu (2014) computed constraints on it from the Cosmic 
Microwave Background (CMB), the large scale structure and Su¬ 
pernovae Type la, while Calabrese et al. (201 ) forecasted errors 
from the CMB on the parameters of an early dark energy possess¬ 
ing anisotropic stress. It is also complementary to that of Amendola 
et al. (2014); Cardona et al. (2014); Sawicki et al. (2013), who also 
put constraints on anisotropic dark energy, but used different mod¬ 
els for its evolution. 

Currently, there are no significant deviations from dark energy 
being a perfect fluid, for instance the result by (Bean & Dore 2004) 
who find c 2 < 0.04 at low significance from CMB-data, such that 
tests whether dark energy is an ideal fluid will be the domain of fu¬ 
ture experiments: Quite generally, the sensitivity to non-ideal cos¬ 
mic fluids requires their respective density to be large enough and 
their equation of state not to be too negative for dark energy per¬ 
turbations to be sufficiently strong (Erickson et al. 2002; Koivisto 
& Mota 2006; de Putter et al. 2010; Ballesteros & Lesgourgues 
2010; Archidiacono et al. 2014). At a first sight it would appear 
that choosing a dark energy equation of state too far from the cos¬ 
mological constant value is incompatible with present constraints 
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(Ade et al. 2015). However, when including extra parameters such 
as the speed of sound and viscosities in the fluid, constraints be¬ 
come much more loose (Mota et al. 2007; Archidiacono et al. 2014, 
See e.g.). 

This article is structured as follows: We develop the nec¬ 
essary perturbation equations for non-ideal dark energy fluids 
and a suitable parametrization in Sect. 2 and discuss cosmo¬ 
logical probes in Sect. 3, before computing forecasts on non¬ 
ideal dark energy properties in Sect. 4. We summarise our re¬ 
sults in Sect. 5. The reference cosmological model is a spatially 
flat, dark-energy dominated model with the parameter choices 
{Q m h 2 , a b lr, n s , D,„) = (0.142, 0.022, 0.67, 0.96, 0.32). This cor¬ 
responds to the constraints from Planck (Ade et al. 2014b) and 
WMAP polarization low-multipole likelihood (Bennett et al. 2013; 
Ade et al. 2014a), and it represents the present official baseline for 
Euclid forecasts. The amplitude of the primordial power spectrum 
was fixed to A s = 2.1 x 10 -9 . The dark energy equation of state pa¬ 
rameter was set to be constant with a numerical value of w = -0.8. 


2 COSMOLOGY WITH NON-IDEAL FLUIDS 
2.1 Expansion dynamics 


Since we focus on late cosmological times, where dark matter and 
dark energy are dominating the energy density of the Universe, we 
can approximate the Hubble function H(a ) = d/a with 

H 2 = H 2 [o m , 0 fl 3 + (1 - n m ,oK 3(1+ “' ) ], (1) 

where a is the scale factor, O m 0 is the dark matter density parameter 
today, w is the equation of state of dark energy, which we assume to 
be constant, and Ho is the Hubble parameter today. In addition, we 
do not consider global curvature. The comoving distance is defined 
as 


x -cf 

a 


da 

a 2 H(a) 



da 

a 2 H(a)/H 0 ’ 


( 2 ) 


with the Hubble distance Xh = c/H 0 ^ 2996.9 Mpc/h. At the 
same time, this defines conformal time r through x = ct. In the 
following, we will set c = 1. 


2.2 Perturbations and their analytical solutions 


If we consider a non-ideal fluid dark energy, characterised by a con¬ 
stant equation of state w, a speed of sound c s , and an anisotropic 
stress component <x, we can write the evolution of cr as in Hu 

(1998): 

,38 cl V 

a + ~a a ~ 3 (1 +w) 2 ~cP r H' ( 


where the prime indicates derivative with respect to a and c 2 is 
called viscosity parameter, as it gives a measure of the fluid’s 
viscosity: Indeed, Eq. (3) implies that when c 2 = 0, then the 
anisotropic stress component cr is also vanishing, while when e.g. 
cj = 1/3 the evolution of anisotropic stress for radiation up to the 
quadrupole is recovered. 

To this equation, we add the first order perturbation equations 
for the density contrast 6 and the velocity perturbation V 


6 ' 

V' 
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where 5p is the pressure perturbation, p is the dark energy den¬ 
sity, 0 and 0 are the metric perturbations in the Newtonian gauge, 
defined by the line element 


dv 2 = a 1 [—(1 + 20)dr 2 + (1 - 20)d.r,dv'J. 
Pressure perturbations are parametrized as 


, 3aH(c 2 - c 2 ) 

Sp = c\pS + -^- ^pV, 


( 6 ) 


(V) 


where c 2 = pip = w is the adiabatic speed of sound for a fluid 
with constant equation of state, to which c 2 reduces in the case 
of a perfect fluid, when no dissipative effects, leading to entropic 
perturbations, are present (Bean & Dote 2004). 

In order to close the differential equation system, one needs to 
include the Poisson equation 


^0 = -47rGfl 2 ^ pi |<5, + ^- Vi 


-4nGa 2 ^ p, A,, 


( 8 ) 


(where the sum runs over all clustering fluids, G is the Newton con¬ 
stant, and in the last equality we have defined the gauge-invariant 
density perturbation of the i-th fluid, A, = <5 ; + 3aHVi/k 2 ) and the 
fourth Einstein equation 


\2jiGa 2 (1 + w)pcr 

(9) 

9 -Hl{\ - G m , 0 )a- (1+3w >(l + w)cr 


B(a) tr. 
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In Sapone & Majerotto (2012) the following analytical solutions 
for 8, V and tr were found for the matter dominated era: 
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(13) 


where k 2 0 = -0o, which is valid strictly only during matter domi¬ 
nation and while neglecting dark energy perturbations. 

As found in Sapone & Majerotto (2012), to which we refer for 
further detail on the analytic solutions, the relevant quantity is the 
effective sound speed 
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2 

eft 
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(14) 


as Eqs. (11-13) can be rewritten in terms of it. This means that 
the sound speed and the viscosity have a similar damping effect 
on density and velocity perturbations (as also noticed in Mota et al. 
2007; Calabrese et al. 201 ). It is interesting to notice that the effect 
of c 2 is enhanced with respect to that of c 2 by a factor of 8(c 2 - 
w)/[3(l + w)], which is ~ 10 if w ~ -0.8 and for very small c 2 , 
which are the cases where a viscosity can be observed best, as will 
be shown in the following sections, and bounded by the case of a 
cosmological constant, as Eq. (11) diverges for w = — 1. 


2.3 Observable parameters 


To understand how the viscosity affects the physical observables, it 
is useful to introduce the clustering parameter Q and the anisotropy 
parameter p, defined in Amendola et al (2008) and computed in 
the case of viscous dark energy in Sapone & Majerotto (2012). 

Q parametrizes the deviation from a purely matter-dominated 
Newtonian potential and is given by (see Sapone & Majerotto 2012) 


Q -1 = 


Pm A,„ 


i - a 


•m, 0 




(1 + w) 


1 -3w + -J^-C 2 




= Qo 


1 + a a’ 


(15) 


where a = 2fc 2 c 2 ff /[(3//5f2,„, 0 )(l - 3w)] and go = (1 + w)(l 
f^m.o)/ [f2m,o(l - 3vt’)]. 

The anisotropy parameter is then given by 


r,= ^-l=- -H 2 (l-n„ 1 > 0 )(l+w)- 


k 2 Q 



(16) 


This is nonzero only when anisotropic stress is present and the met¬ 
ric perturbations 0 and 0 are different. 

Let us finally define the parameter Z (Amendola et al. 2008) 
as 


2 = 11 + -Z7| <2. 


(17) 


This is useful because it represents the deviation of the weak lens- 
ing potential ® = 0 + 0 from its behaviour in the case of no dark 
energy perturbations. 


3 COSMOLOGICAL PROBES 

In Sapone et al. (2013), we forecasted constraints to the viscosity 
parameter and the sound speed from the Euclid galaxy clustering 
and weak lensing surveys. Here, we aim to complete the picture by 
adding to the latter the constraints from the iSW tomography signal 
obtained by cross-correlating galaxies mapped by the Euclid photo¬ 
metric instrument with the Planck temperature map: This provides 
a combination of all major probes of cosmic structure formation, 
which draw their sensitivity from the growth rate and interactions 
between fluids, from the shape of the initial perturbations and from 
the expansion history. 

The Euclid survey is a mission of the ESA Cosmic Vision 
program that will be launched in 2020, and will perform both a 
photometric and a spectroscopic survey, the first aiming mainly at 
measuring weak lensing while the second at measuring the galaxy 
power spectrum. The Planck satellite is also a mission of ESA Cos¬ 
mic Vision program, already operating and mapping the CMB fluc¬ 
tuations with unprecedented precision and control of systematic ef¬ 
fects. 

To perform our forecasts we use the Fisher matrix (Tegmark 
et al. 1997), which quantifies the decrease in likelihood if a model 
parameter 9 a moves away from the fiducial value, and can be 
computed for a local Gaussian approximation to the likelihood 
JL on exp(-^ z /2). In our forecasts we assume the official Euclid 
specifications, that can be found in Laureijs et al. (20 II). The fidu¬ 
cial cosmological parameters correspond to the 2013 Planck mea¬ 
surements (Ade et al. 2014b), except for the value of w, for which 
we assume w = -0.8, in order for the effects to be more clearly 
visible (as done in Sapone & Majerotto 2012; Sapone et al. 2013) 
and of course except for the values of c 2 and c 2 . 
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In the following, we will describe the iSW tomography signal 
and give a short summary on the signal coming from weak lensing 
and the galaxy power spectrum from spectroscopy. 


3.1 iSW signal 


When a CMB photon moves into a time-evolving metric such as 
that of Eq. (6) the unbalance between the blue-shift experienced at 
the entrance and the red-shift experienced at the exit of its varying 
potential well originates a perturbation £ in the CMB temperature 
?cmb given by (Sachs & Wolfe 1967), 


( = 


AT 

Tcmb 



dc/> di//\ r XH 

+ Jo 


A\a 2 H 


<9® 
da ’ 


(18) 


where x is the comoving distance (see Eq. 2) and ® is the weak 
lensing potential. 

In the case of pure matter domination, ® = const, hence the 
iSW effect vanishes, while in presence of any fluid with w 0 the 
temperature fluctuation £ will be nonzero, so that the late iSW is 
particularly interesting to us, as it is originated by the appearing 
of dark energy and it is an independent proof of its existence (first 
detected by Boughn & Crittenden 2004; Giannantonio et al. 2012, 
for an updated measurement). 

Let us now compute the term inside the integral, passing to 
Fourier space, and in the case of viscous dark energy (see also 

Sapone & Majerotto 2012; Schaefer 2009): 


<9® 

da 


3 H- 0 n m , 0 

2 ak 2 


£ (a, k) A' (a, k) + 


+ £' ( a , k) A,„ (a, k) -£ (a, k ) A,„ ( a , k) >. 

a 


(19) 


It is possible to see from this expression that anisotropic pertur¬ 
bations enter the iSW effect in two ways: by modifying A„, and 
through the additional presence of £ and £'. At linear order, it is 
possible to isolate today’s A,„ from its time evolution: 


A,„ (a, k) = aG(a, k) A m 0 ( k ), 


( 20 ) 


where A m o ( k ) = A m (a = 1, k). We write hence Eq. (19) as: 


3® 3 Hlkl m0 <9 ( i 

7ki = ~2 da r {a ’ k) E (a ’ k) ) Ara -° {k) ’ 

and Eq. (18) reads now 

r><H 

C = A x W ( (x)A m , 0 {k) 

Jo 

where the weighting function Wfy) I s 


W ( Ct) = 


3 o 2 d 

r Cl tl ~~ 

2 kr da 


|G(a,/fc)£(a,k)|. 


( 21 ) 


( 22 ) 


(23) 


Since the iSW is a secondary effect of the CMB (Rees & Sciama 
1968), it can be separated through cross-correlation to the galaxy 
density (Crittenden & Turok 1996). Let us hence write the galaxy 
density obtained through imaging surveys, in order to compute its 
cross-correlation with the iSW. The line-of-sight projected galaxy 
density y is given by (Smail et al. 1995) 


7= [ X " & X D(z)^b(x)G(x)8(z) (24) 

Jo ax 

being D(z) the galaxy distribution defined as 


D(z) = 



exp 



(25) 


where the Hducial parameters p D and zo depend on the imaging 
survey considered. In the case of Euclid, they are p D = 3/2, zo = 
ttmean/ V2, and z mean = 0.9 (Laureijs et al. 201 1). 

Even though the signal from the iSW increases noticeably 
when cross-correlating it with the galaxy density field, both the 
cross-correlation spectrum and the galaxy spectrum are line-of- 
sight integrated quantities, hence much information may be lost. 
For this reason we decide to use iSW tomography (Ho et al. 2008; 
Douspis et al. 2008; Juergens & Schaefer 2012), and in particular 
we divide the whole galaxy sample into 5 bins with equal num¬ 
ber of galaxies (in order to match with the binning used by official 
Euclid documents for weak lensing tomography). To do this, we re¬ 
place the galaxy distribution function D(z), Eq. (25) in y, Eq. (24), 
with the radial distribution function of galaxies in the z-th bin D,(z), 
obtained by binning the overall distribution D(z) and convolving it 
with the photometric redshift distribution function (Amendola et al. 
2008). 

We are finally able to write our observable, i.e. the iSW-galaxy 
cross-correlation spectrum Cf r ,;(T) in the z-th redshift bin, along 
with the iSW-auto correlation spectrum C^{C) and the galaxy- 
galaxy auto correlation spectrum C yyi j(£) of the ij- bins, (which 
are both needed in order to estimate statistical errors coming from 
C(((£)), by applying a Limber-projection (Limber 1954) in the flat- 
sky approximation: 


II 

f * % W ( (x)W y ,, ; (x) P AA (k = t! X \ 

Jo X 

(26) 

C a (€) = 

r x H dy 

4 wl or) p aa (k = e/ X ) 

(27) 


Jo X~ 


Cyy.,j(0 = 

f " % w y ,i(x)w yJ (x) p AA (k = e/x) 

(28) 


where P AA Ik) is the linear matter power spectrum today, and the 
galaxy weighting function of the z-th bin W y j( X ) is 

W yJ (x) = D,(z)^b(x)G(x). (29) 

The tomographic iSW spectra Q y- ,(T) are shown in Fig. 1 for two 
fiducial models: a standard dark energy model with cl = 1 and 
cl = 0 (solid lines) and a model with viscosity: cl = 1CT 5 and 
cl = 1(L 6 (dashed lines), where the colour shading indicates the 
reshift bin for which C( y j(£) was evaluated. The iSW-effect is a 
large-scale effect originating from low redshift, as the influence of 
dark energy on the growth of gravitational potentials in the large- 
scale structure is strongest. The effect of dark energy visosity and 
small sound speed is strongest on large scales as well (see also 
Sapone & Majerotto 2012), and affects a wide range of multipoles. 
Keeping all cosmological parameters fixed, dark energy viscosity 
would increase the amplitude of the iSW-effect by up to 25% on 
large angular scales and at low redshift. This sensitivity of the spec¬ 
tra at low multipoles is fortunate because these scales can be well 
probed with the iSW-effect. 


3.2 Weak Lensing 

To the iSW tomography signal we add the weak lensing tomo¬ 
graphic signal (Hu 1999,2002; Heavens 2003; Jain & Taylor 2003), 

coming from the same photometric survey as y and using the same 
redshift bins. Here, we only give the main equation expressing the 
weak lensing power spectrum, which is used for our forecasts, and 
refer to Sapone et al. (2013) for further details. 

In presence of anisotropic stress, the weak lensing conver- 
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Figure 1. Tomographic iSW-spectra Cr 7 ,i(C) for two dark energy models: 
c 2 = 1 and c 2 = 0 as well as c 2 = 10 , c 2 = 10 6 - Blue to light green lines 
correspond to redshifts z in the intervals [0.01 - 0.5595] (blue). [0.5595 - 
0.7871], [0.7871 - 1.0165], [1.0165 - 1.3184] and finally [1.3184 - 2.5] 
(green). 


Figure 2. Tomographic weak lensing spectra C KK ji{() for two dark energy 
models, c s = 1 and c v = 1 as well as c s = 10 -5 , c„ = 10 ( \ both includ¬ 
ing the shape noise term. Blue to light green lines correspond to the same 
redshift binning as in Fig. 1: [0.01 - 0.5595] (blue), [0.5595 - 0.7871], 
[0.7871 - 1.0165], [1.0165 - 1.3184] and finally [1.3184 - 2.5] (green). 


gence power spectrum is given by (Jain & Taylor 2003; Hu 2002, 
1999; Hu & Jain 2004) 

d y 

4 W Kj (x)W Kj (x) £ 2 P N L(k = £/x,X)- ( 30 ) 

o X 

where the subscript ij refers to the redshift bins around Zi and Zj, 


with 

3,0, m 

Fi(x ) 


W Kj (x) = 

24 

„ X 
a 

(31) 

PM = 

nxn 

dx n(x ) A Cl ) 

X 

(32) 


and where Dj is the same tomographic distribution function of 
galaxies used for the iSW-effect. While tomography in general 
greatly reduces statistical errors the actual shape of the choice of 
the binning does not affect results in a serious way, although in prin¬ 
ciple there is room for optimisation (Schafer & Heisenberg 2012). 

In Fig. 2 we show the tomographic weak lensing spectra 
C K ji(t) for the same models and the same redshifts as in Fig. 1. 
As for iSW, the effect of viscosity is detected at large scales and 
for a large range of scales (but smaller than for iSW). Instead, con- 
trarily to iSW, here the sensitivity to viscosity is stronger at higher 
redshift. This is because the efficiency of weak lensing is higher for 
longer light paths. 

In principle it is also possible to define a cross-spectrum of 
weak lensing and iSW, and of weak lensing and galaxy 

distribution, C^yj(P), but both these spectra are subdominant with 
respect to C( y {£). This is because the weak lensing convergence 
signal comes from the distortion of the light path at redshifts in¬ 
termediate between us and the galaxies mapped by the imaging 
survey, while the iSW signal originates precisely at the same red¬ 
shifts where the galaxies are. We have tested this fact by computing 
the signal to noise-ratio for measuring Cz K j{k) and found it much 
smaller than that of Qy.^C). 

3.3 Spectroscopic galaxy power spectrum 

To the iSW and weak lensing measurements, both measured 
through photometric observations, we add data coming from the 


power spectrum of spectroscopically observed galaxies. Here we 
only show the expression of the observed power spectrum, which 
is needed in order to compute our forecasts, and refer again the 
reader to Sapone et a (20 1 3) for further detail. 

Following Seo & Eisenstein (2003) we write the observed 
galaxy power spectrum as: 

P^(.Z,k r ,Ur) = G2 (Z,k)b(z) 2 (l + Pfl 2 ) 2 Por(k)+P shot , (33) 

where the subscript r refers to the reference (or fiducial) cosmolog¬ 
ical model. 

Here P shot is a scale-independent offset due to imperfect re¬ 
moval of shot-noise, /j = k ■ r/k, is the cosine of the angle of the 
wave mode with respect to the line of sight pointing into the di¬ 
rection r, P 0r is the fiducial matter power spectrum evaluated at 
redshift zero, G(z, k) is the linear growth factor of the matter per¬ 
turbations, b(z) is the bias factor and Da(z ) is the angular diam¬ 
eter distance. The wavenumber k and ^ have also to be written in 
terms of the fiducial cosmology (see Seo & Eisenstein 2003; Amen- 
dola et al. 2005; Sapone & Amendola 2007, for more details). The 
fiducial bias used can be found in Orsi et al. (2010), who derived 
their results by using a semi-analytical model of galaxy formation, 
while the matter power spectrum has been computed with a modi¬ 
fied version of the CAMB code 2 (Lewis et al. 2000) accounting for 
anisotropies. 


4 STATISTICAL ERRORS FORECASTS 

In this section we estimate marginalised statistical errors on the 
sound speed and viscosity parameters c 2 and c 2 through the Fisher- 
matrix formalism (Tegmark et al. 1997), which assumes a Gaussian 
likelihood and unbiased measurements. 


2 http://camb.info 
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4.1 iSW Fisher matrix 


4.3 Spectroscopic galaxy distribution Fisher matrix 


The sensitivity of line of sight-integrating effects can be boosted by 
subdividing the galaxy population into redshift bins: For the iSW- 
effect this was first carried out successfully by Ho et al. (2008), and 
systematically investigated by Juergens & Schaefer (2012). 

The Fisher-matrix of the iSW-effect follows directly from the 
variance of the spectrum estimates, 


E-iSW 
b afi 


l 


dC c yi{t) , 

^g-Cov-H) 


dC (y ,j{t) 

dBp 


(34) 


where the sum runs from C = 5 to £ = 300 3 , 8 a are the cosmological 
parameters, CoVy(£) is the covariance of the spectrum C^ y j(i) and 
is given by 


The galaxy power spectrum Fisher matrix is given by (Seo & Eisen- 
stein 2003) 

gc f W Edk dinP^ c (z; k,n) d InP^ ec (z; k ,//) 

F * = J^ 4 ^ dB a W, XVeff ’ (41) 

where GC stays for galaxy clustering, the observed galaxy power 
spectrum is given by Eq. (33), the derivatives are evaluated at 
the parameter values of the fiducial model, k min = 0.001 and k m . dX is 
such that the rms amplitude of the fluctuations at the corresponding 
scale i? max = 2 n/k mdx is <x 2 (R max ) = 0.25, with an additional cut at 
kmax = 0.2/7/Mpc. in order to remain in the linear regime, V eff is 
the effective volume of the survey, given by 


Covy (£) = — l — f C {7 jC {7 j(t) + C a (nc yy jj({) 1, (35) 

+ 1 /sky 

and where quantities with the bar represent the estimate of the sig¬ 
nal, including intrinsic CMB fluctuations, instrumental noise and 
the beam of the CMB experiment as noise sources: 

C*w(0 = C Cy ,m (36) 

Cfftf) = Cg(£) + CcmbW + wf B~ 2 (€) (37) 

C yyJj (0 = C nij (f)+ ^ (38) 


For Planck’s noise levels, Wj, 1 = (0.02/rK) 2 has been used 
and the beam was assumed to be Gaussian, B~ 2 (P) = 
(2 x 10 -8 ) 2 exp[A0 2 £(i + 1)], with FWHM-width of A 9 = 7'.1, 
corresponding to channels of Planck closest to the CMB-maximum 
at ~ 160 GHz. n, is the number of galaxies per steradian in the 
tomography bin i. We assume uncorrelated noise terms, and as a 
consequence the cross-spectra C( y j{£) are unbiased estimates of the 
actual spectra, see Eq. (36). The spectrum Ccmb( 0 of the CMB pri¬ 
mary anisotropies from Planck has been computed with the CAMB 
code. 


4.2 Weak Lensing Fisher matrix 


The Fisher matrix for weak lensing is given by: 


a(i 


' /sky ^ 


(2f + 1) dC KKjj ({) dC KK , km {t) , 


dB a 


Cf(E- 


dB B 


-cftf) 


(39) 


where the sum runs from t = 5 to £ = 5000, (as from the offi¬ 
cial Euclid prescriptions, see Laureijs et al. 201 1), and where the 
sum over repeated indices is implied. We added a Poissonian shape 
noise term to the weak lensing spectra, 


C KK ,ij(t) - C KK jj(£) + Sjj 



(40) 


■yin, is the rms intrinsic shear (here, we assume {y l Jp)= 0.22) and «,• is 
the number of galaxies per steradians belonging to the i-th bin. We 
assume a Gaussian shape of the covariance while noting that non- 
Gaussian contribution can have a strong influence on the derived 
forecasts (Takada & Jain 2009). 


3 The integration range for the iSW-effect as well as the details of instru¬ 
mental noise and angular resolution are not very important, as most of the 
signal is at low t below t ~ 100, due to the large cosmic variance provided 
by the primary CMB fluctuations, which is the largest source of noise. 


/ nPffiEp) \ 

U/>r(M)+ij y 


2 


(42) 


the latter equation holding for an average comoving number density 
Pi. The number densities and further fiducial Euclid specifications 
can be found in Laureijs et al. (2011), Majerotto et al. (2012). 


4.4 Forecasts 

We computed forecasts on the measurement of cj and c 2 for a wide 
range of fiducial values, in order to capture the parameter determin¬ 
ing capability of both experiments for a previously unknown set of 
parameters. The probes are assumed to be uncorrelated as discussed 
above, hence their Fisher-matrices add, 

~ b aP + b a p • I 43 ) 

and we derive confidence contours on c 2 and c 2 . and individual er¬ 
rors from this combined Fisher-matrix, marginalising over all other 
five parameters considered in this analysis. has been further 
marginalised over P shot , while the galaxy bias has been kept fixed. 

Our forecasts on the following fiducial models: {c 2 ,c 2 } = 
jl, 0), {10 3 , 1(T 4 ), {10 5 , 10 6 } and {10 6 , 1(T 6 ) are shown in 
Figs. 3, 4, 5 and 6, respectively. The first model corresponds to 
the case of simple scalar field dark energy, while the following two 
pairs of fiducial models were chosen such that c 2 = 10 c 2 because, 
as mentioned previously in Sec. 2.2, the relevant quantity is the ef¬ 
fective sound speed, and the effect of c 2 in it is — 10 times stronger 
than that of c 2 when w = -0.8 because of the factor multiplying c 2 
in Eq. (14). The last model, also having small c 2 and c 2 , does not 
verify the latter relation, and has been chosen in order to be com¬ 
pared to previous work (Sapone et al. 2013) and to the similar case 
{c 2 ,c 2 ) = {10- 5 , 10 6 }. 

In all plots, iSW constraints are shown in blue, weak lensing 
ones in dark blue, GC ones in green, combined iSW-GC ones in 
yellow, combined iSW-weak lensing ones in orange, and combined 
iSW-GC-weak lensing ones in red. 

From Figs. 3-6 it is clear that the results depend very much 
on the chosen fiducial model. A common feature is that iSW on its 
own does not provide very strong constraints. In particular (see also 
Tab. 1), c1 is quite badly constrained, with relative errors ranging 
between 6.1 x 10 4 and 1.1 x 10 5 , while relative errors on c 2 are 
much smaller: between 1.1 x KT 1 and 1.4 x 10. This was to be ex¬ 
pected since the iSW-effect along has a rather small signal strength 
of about 5cr for cross-correlating the CMB with the Euclid galaxy 
sample (Douspis et al. 2008). 
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Figure 3. Forecasted lcr-constraints on c 2 and c 2 for individual probes and 
all possible combinations, for the fiducial choice c 2 = 1 and c 2 = 0. 


Figure 4. Forecasted lcr-constraints on c 2 and c 2 for individual probes and 
all possible combinations, for the fiducial choice c 2 = 10 3 and c 2 = 10 4 . 
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1 

0 
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WL 

10- 3 
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1.4 x 10 2 


10- 5 
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1 

0 

1.1 x 10 2 

cr.2 = 9.1 

GC 

10- 3 

10- 4 

7.5 x 10~ 2 

7.2 X 10~ 2 


10- 5 

10- 6 

2.2 

1.9 


10- 6 

10- 6 

4.7 

1.2 


1 

0 

2.9 x 10- 1 

<r c 2 = 2.7 

iSW+WL 

10- 3 

10- 4 

2.5 x 10- 1 

5.6 


10" 5 
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2.6 x Hr 2 

7.5 x 10 1 


10~ 6 

10- 6 

7.0 

7.6 x 10~‘ 


1 

0 

7.2 x 10~ 2 

cr c 2 =3.0 

iSW+GC 

10- 3 

10- 4 

5.3 x 10- 2 

5.5 x 10~ 2 


10" 5 

10- 6 

1.6 x 10” 2 

1.2 
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3.9 

4.2 x 10 1 
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0 
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<r c 2 =3.7 

GC+WL 
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7.4 x 10~ 2 

7.1 x 10~ 2 
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1.9 

1.7 


10~ 6 

10- 6 

3.5 

4.8 x 10~‘ 


1 

0 

6.7 x 10~ 2 

cr c 2 = 2.4 

all 

10- 3 

10- 4 

4.5 x 10~ 2 

4.5 # x 10~ 2 

10- 5 

10- s 

1.2 x 10~ 2 

6.1 x 10- 1 


IQ- 6 

10- s 

3.3 

3.4 x 10~‘ 


Table 1. Relative errors on the parameters c~ and c 2 from iSW, weak lensing 
and GC alone, from the combination of iSW and WL, iSW and GC, GC and 
WL, and from all three datasets. For the case cl = 0 the absolute error cr 2 
is given. 


Weak lensing constraints 4 are much stronger than iSW ones 
in the case of a fiducial scalar field dark energy, but become pro¬ 
gressively comparable to them when the fiducial cl and c 2 become 
smaller, with the exception of the case c 2 = c 2 = KT 6 . 

Even though both iSW and weak lensing do not give very 
strong constraints on sound speed and viscosity (see also Tab 1), 
it is very interesting to notice that the two data sets complement 
each other very well. This is especially true for the case cj = 10~ 5 
and cl = 10 -6 , represented in Fig. 5, where the blue ellipses, which 
indicate errors from iSW, have a very different degeneration direc¬ 
tion with respect to the dark blue contours, corresponding to errors 
from weak lensing, but are comparable to them in size. Therefore 
the resulting combined errors are much smaller than those from a 
single dataset. In particular, the iSW effect gives better constraints 
on the sound speed and weak lensing on the viscosity parameter. 

Also, in the case of Figs. 3, 4 and 6 iSW and weak lensing have 
different degeneracies, but here joining them does not improve the 
errors significantly because the weak lensing effect gives stronger 
constraints on both parameters. Here, the improvement in combin¬ 
ing the two probes is rather the multiplication of a constraining 
likelihood with a wide one, resulting nevertheless in an increase in 
peakiness. 

Also errors on cl and cj from GC are orthogonal to those from 
the iSW. but only in the case where cl = 10~ 5 and cl = 10 G , see 
Fig. 5, this helps reducing the errors, because the former dataset 
performs better in constraining cj and the second cl. In the case of 
fiducial cl = 10 -3 and c; = 10~ 4 it is GC which gives best errors 
on both parameters, while for the cl - 1 and cl = 0 fiducial model, 
it is weak lensing. 

Another interesting question is whether iSW adds important 
information to that provided by the other two datasets, which had 
already been analysed in Sapone et al. (2013). Table 1 answers this 
question. It turns out that the information from iSW helps signif¬ 
icantly in constraining cl and cj if the true model has cj = 10~ 5 
and cl = 10 -6 , see Fig. 5. In this case the iSW alone gives a strong 
constraint, which has moreover a different degeneration direction 
with respect to the error from galaxy clustering and weak lensing. 


4 With respect to Sapone et al. (2013) we have improved the estimation of 
Pnl by using the full CAMB output instead of an analytical approximation 
to it. 
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Figure 5. Forecasted lcr-constraints on c ^ and cl for individual probes and 
all possible combinations, for the fiducial choice cl = 10 -5 and cl = 10 -6 . 
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Figure 6. Forecasted lcr-constraints on c 2 and c 2 for individual probes and 
all possible combinations, for the fiducial choice c 2 = 10 6 and c 2 = 10 . 


For the other three fiducial models the gain when adding iSW is not 
very strong, as in both cases the combination of weak lensing and 
galaxy clustering gives a much tighter constraint in both the sound 
speed and the viscosity than the iSW alone. 

It is interesting to notice (see Fig. 6) that when both the sound 
speed and the viscosity are small, but the relation between c 2 and c 2 
differs from c 2 ~ lOcJ, the iSW effect error ellipse becomes much 
larger and as a result the sound speed parameter is less strongly con¬ 
strained. Thus we conclude that (i) very interesting results can be 
obtained through a combination from different cosmological probes 
and that the iSW-effect is able to tighten constraints signifi¬ 
cantly for cases where there is a large difference between c 2 and 


5 SUMMARY 

In this paper we have investigated how well the viscosity and 
sound speed of dark energy can be measured with the iSW cross¬ 
correlation spectrum, when using Planck and Euclid observations. 


and how joining iSW measurements to galaxy clustering and weak 
lensing ones improves constraints. 

We found that the speed of sound is quite well constrained, 
with relative errors as small as 0.14 for small fiducial c 2 and cj, 
while relative errors on the viscosity parameter are very large. Even 
though the anisotropic stress is not well constrained by the iSW, the 
error ellipses are interestingly orthogonal to those from weak lens¬ 
ing, hence the combination of these two datasets constrains tightly 
the parameter space, giving relative errors on c 2 and c 2 as small 
as 2.6 x 10~ 2 and 7.5 x 10~ 2 respectively. This is an improvement 
of a factor ~ 1500 in the measurement of the sound speed and 
~ 50 in the measurement of the viscosity parameter, with respect to 
the weak lensing only constraint. The improvement obtained when 
combining iSW with galaxy clustering is smaller: a factor of ~ 1.5 
in cl and ~ 150 in c 2 . Finally, the addition of iSW to weak lensing 
and galaxy clustering constraints is most important if the fiducial 
sound speed and viscosity parameter are very small, while it is not 
very relevant for higher fiducial values of c 2 . 

It is also important to remind that in order to make the ef¬ 
fect of dark energy perturbations stronger we have always used 
a value of the equation of state parameter w = -0.8. For val¬ 
ues close to w = -1 the effects on the observables due to the 
dark energy perturbations are reduced, as all the phenomenologi¬ 
cal functions used (such as Q(k,a)) have a term oc (1 + w). If we 
use a value of w = -0.9 we expect our final errors on the param¬ 
eters to increase. But by how much? All the observables used in 
this paper depend most strongly on Q 2 (see Eq. 15) which is in¬ 
trinsically included into the matter power spectrum; for a sound 
speed equal to zero Q — 1 = (1 + w)/(l - 3w)a~ 3w so the rel¬ 
ative increase of the errors on the sound speed will be given by 
1 /[(<2(w = -0.9) - 1 )/(Q(w = -0.8) - l)] 2 which is of about a 
factor 4 larger, in agreement also with the results found in Sapone 
etal. (2010). 

A detection of sound speed and viscosity different from the 
values associated to a classical scalar field, i.e. c 2 = 1 and cj = 0, 
will point to a new understanding of the accelerated phase of the 
Universe. This is because the non ideal fluid considered in this pa¬ 
per can be thought of as en effective dark energy fluid parametriz¬ 
ing a modified gravity model, see Kunz & Sapone (2007). In prac¬ 
tice, the detection of a zero sound speed does not automatically 
mean that we are dealing with an actual dark energy fluid, even 
though one would nevertheless experience effects which could be 
attributed to fluctuations of a fluid. 

In this paper we found that joining data from Euclid and 
Planck we are able to constrain simultaneously the sound speed 
and the viscosity parameters, provided that the two are sufficiently 
small. This is mostly due to the different sensitivity of the three 
observables, i.e. GC, WL and iSW to the two parameters. In most 
cases the iSW has a different degeneracy with respect to WL and 
GC, and this helps reducing the errors on c 2 and c 2 by a factor of 
~ 100 (as pointed out before). Our results are in agreement with 
what found by Mota et al. (2007); Calabrese et al. (2011); Chang 
& Xu (2014), who show that for values of c 2 approaching c 2 = 1 
the detection of a positive viscosity is very difficult, even when, as 
in the case of Calabrese et al. (201 1 ), an early dark energy helps its 
detection by increasing its effect on smaller scales. 

To conclude, are Euclid and Planck able to measure the sound 
speed and the viscosity parameters of a dark energy component? If 
the values of c 2 and cj are small enough, the answer is yes; conse¬ 
quently, we will be able to constrain well the effective dark energy 
model. On the contrary, if sound speed and viscosity will escape 
detection, at least one of the two parameters will likely have large 
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values. We would assume that other cosmological probes would not 
directly provide constraints on dark energy properties, but would 
nevertheless be able to provide constraining power by fixing other 
parts of the cosmological model, such as the dark matter density or 
the dark energy equation of state, which was not subject to variation 
in our investigation. 
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